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Lattice simulations on SU(2) and SU(3) gauge theories with matter fields in the fundamental, 
adjoint and two index symmetric representations are needed to determine if these theories are near 
or within the conformal window as required for their applications in beyond standard model phe- 
nomenology. Simulations with Wilson fermion action are subject to artifacts linear in the lattice 
spacing a, and must be improved. We provide the necessary coefficients for perturbative improve- 
ment of the boundary terms when using Schrodinger functional boundary conditions and furthermore 
show that correctly implemented O(a) improved actions are necessary for reliable results. 
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There has recently been interest in studying quan- 
tum field theories with a nontrivial infrared fixed point 
(IRFP), both in continuum and on the lattice. Un- 
der the renormalization group evolution, the coupling of 
these theories shows asymptotic freedom at small dis- 
tances, analogously to QCD, but flows to a fixed point 
at large distances where the theory hence looks confor- 
mal. Such theories have applications in beyond Standard 
Model model building. These include unparticles, i.e. an 
infrared conformal sector coupled weakly to the Standard 
Model and (extended) technicolor scenarios, that ex- 
plain the masses of the Standard Model gauge bosons 
and fermions via strong coupling gauge theory dynamics 
[343] - In addition to direct applications to particle phe- 
nomenology, the phase diagrams of gauge theories, as a 
function of the number of colours, N, flavours Nf and 
fermion representations, are interesting from the purely 
theoretical viewpoint of understanding the nonperturba- 
tive gauge theory dynamics from first principles. 

Several methods to estimate the vacuum phase dia- 
gram of a gauge theory exist. A traditional method is 
the ladder approximation to the Schwingcr-Dyson equa- 
tion for the fermion propagator yielding an estimate for 
the onset of chiral symmetry breaking and signaling the 
departure from conformal to confining phase [3, 0] . How- 
ever, the only truly first principle method is constituted 
by lattice simulations. Several initial studies have ap- 
peared in literature: for example SU(2) with fundamen- 
tal representation fermions 0, Hj], SU(2) with adjoint 
fermions | 9rjl4ll and SU(3) with fermions in the funda- 
mental (l5Hl?| or in the two-index symmetric (l8j . i.e. 
the sextet, representation. 

The studies with Wilson fermions are subject to lat- 
tice artifacts proportional to the lattice spacing a. In 
the context of high precision results for QCD like the- 
ories, i.e. SU(N) gauge theory with modest number of 



flavors, a program to cancel these lattice artifacts has 
been devised 3, 2(j. The basic idea in this approach is 
to introduce local counterterms whose coefficients will be 
fine tuned to cancel all O(a) contributions. In this Let- 
ter we provide the values of the perturbative coefficients 
required in the analyses of SU(2) and SU(3) gauge theo- 
ries with fermions in the fundamental, adjoint and two- 
index symmetric representations. In addition to these 
results, wc highlight their significance: The improved ac- 
tions neglecting some of the counterterms will not work 
and, moreover, the boundary conditions for the back- 
ground field when using Schrodinger functional, must be 
chosen carefully. This point has been emphasised also in 

The basic Wilson lattice action is So = Sg + Sf, where 
Sq denotes the standard Wilson plaquette action and 
5f is the usual Wilson fermion action for Nf (mass de- 
generate) Dirac fermions in the fundamental or higher 
representation of the gauge group. On lattices with pe- 
riodic boundaries, the 0(a) discretization errors in this 
action can be removed by introducing the Sheikholcslami- 
Wohlert -term fTl 



Smxpi = S + a 5 c sw ^2'ili(x)-a l j, v F fJ , v (x)ip(x), (1) 



and tuning the coefficient c sw so that the 0(a) effects in 
the on-shell quantities cancel. Here a M „ = i[7 M ,7„]/2 and 
F fll/ (x) is the symmetrized lattice field strength tensor. 

The Schrodinger functional scheme is often used to 
measure the evolution of the coupling constant. In this 
scheme new contributions to linear lattice artifacts arise 
due to 0{a) errors arising from the fixed spatially con- 
stant boundary conditions at times t = and t = L. 
Here L 4 is the volume of the lattice, and the spatial link 
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variables at the t — and t = L boundaries arc fixed to 



Uk{x)\( Xo=0 ) = cxp(aC fc ), Uk(x)\( Xo =L) = exp(aCfc), 



(2) 



where a is the lattice spacing and 



C fe = — diag((/>i, 



with the constraint 5Zi=i 0« = Si=i $ = 0. For SU(2) 
and SU(3) the boundaries can be parametrized as p2T[23j 

(<h,<h) = (v,-v), (4>i,4>2) = {p-v,-p+v) su(2) (4) 



^ = z diag(^i,...,^ c ), 
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= (v-p,-v/2,-V/2 + p) 
) = {- V -'Sp lV /2 + p lV + 2p) 



SU(3). (5) 



The conventional choice of parameters is (p, rj) = (tt, tt/4) 
for SU(2) and (p,rf) = (tt/3,0) for SU(3). The spa- 
tial boundary conditions are taken to be periodic. The 
fermion fields are set to vanish at the t = and t = L 
boundaries and have twisted periodic boundary condi- 
tions in spatial directions: tp(x + Lei) = exp(iw/5)^(x). 
This improves the condition number of the fermion ma- 
trix 24j. The values of <j>i in ((3]) are restricted within 
the so-called fundamental domain in order to guarantee 
a unique least action solution for the background field 

HQ. 

At the classical level the nontrivial boundary condi- 
tions generate a spatially constant chromoelectric field 
and the derivative of the action with respect to rj can be 
easily calculated; it is proportional to the inverse of the 
bare coupling g" 2 . At the full quantum level the coupling 
can now be defined by 23] 



4=f> 

9 or] i 



'dS cl 
drj 



With this preliminary definition, the running of the cou- 
pling is quantified using the step scaling function: 



E(u,s,L/a) = g 2 (g ,sL/a)\ g 2( go ^ L / a - )=u 
= w+(£ 1 ,o + £ 1 , 1 iV / )u 2 . 



(7) 



The second line gives the formula in perturbation theory 
to one loop order, and the fermion contribution is de- 
noted by Ei.i while £ 1-0 stands for the pure gauge one 
loop contribution. To evaluate the se p erturbative con- 
tributions we use the methods in |24 |25| , and choose 
s = 2. The continuum limit of E^i (£i,o) is given 
by the fermionic (bosonic) contribution to the one loop 
coefficient bo = /3o/(167r 2 ) of the beta function, where 
A = H/3iV c - 4/3T(R)N f . In other words, 



lim E 1 , l /(2&o !l ln2) 

L/a^Q 



1. 



0,1 



(8) 



where 6 ,o = HiV c /(487r 2 ) and 6 ,i = N f T (R) / (12ir 2 ) 
corresponding, respectively, to bosonic and fermionic 
contribution in the beta function. 



The O(a) errors originating from the nontrivial bound- 
ary field can be removed by introducing new terms to the 
action and fine tuning the corresponding coefficients so 
that the 0(a) contributions cancel. Complete analysis of 
the necessary counterterms has been presented in [26[ . In 
the case of an electric background field and after setting 
the fermion fields to zero on the boundary, the required 
counterterms are of the form 

6S G>b = ~(ct-l)Y^Tr[l-U(p t )}, (9) 
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pt 



6S F , b = a 4 (c t - 1) ]T[O t (f) - 6' t (x)} 



(10) 



Here U(pt) denotes plaquettes which touch the t = or 
t = L boundaries. The explicit from of the operators Ot 
and 0[ is not needed in the following; for details, see [26| . 
By tuning the coefficients c sw , q, dt to their proper values 
we can remove all 0{a) errors from our action. Pcrtur- 
batively, the coefficient c sw = 1 + O(gg), and it can be 
determined non-pcrturbativcly in full lattice simulations. 
In the following we will concentrate on the one loop order 
perturbative analysis of Ct and Ct, and use c sw = 1. 

The boundary coefficient Ct has a perturbative expan- 
sion of the form c t = 1 + c^g^ + O{go), and similarly for 
c,f To one loop order the coefficient Ct can be extracted 
from the result of 



27| as cf } 



-0.0135(1)6^ for fun- 
damental fermions, and this generalizes to other fermion 
representations simply by replacing the fundamental rep- 
resentation Casimir operator Cf with Casimir operator 
Cr of the representation R under consideration. The 
results for different gauge groups and fermion represen- 
tations are shown in table HI 



The coefficient 



(6) fermionic parts 



can be split 



tion Cj 1 ' -* is entirely due to gauge fields and has been 
evaluated in for SU(2) and in for SU(3). The 
fermionic contribution c^ 1 ' 1 ^ has been evaluated for fun- 



c^N f 



into gauge and 
The contribu- 



damental fermions in [24j both for SU(2) and SU(3). 
We have extended these computations for SU(2) and 
SU(3) gauge theory with higher representation fermions 
[j^ . The results for the nonzero improvement coeffi- 
cients are tabulated in table HI For the details of the 
numerical method used to determine coefficient c^ 1 , 
we refer to the original literature where the method 
was developed for fundamental representation fermions 
2ol . 24|]. Our results are numerically consistent with the 
generic formula c t (1,1) « 0.019141(2T(i?)), where T(R) 
is the normalization of the representation R, defined as 
Tr(Tgr£) = T(R)S ab . The numerical results collected in 
Table U should be useful for studies of Wilson fermions 
in higher representations, and we highlight two issues: 
First, improvement should be carried out in order to ob- 
tain reasonably accurate results on lattices of reasonable 
size. Second, the improvement program has to be carried 
out consistently and correctly. 
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-0.0543(5) 


0.0192(2) 


-0.0101(3) 
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-0.0543(5) 
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-0.0270(2) 
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3 


-0.08900(5) 
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-0.0180(1) 
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-0.08900(5) 


0.1148(3) 


-0.0405(3) 


3 
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-0.08900(5) 


0.09571(2) 


-0.0450(3) 



TABLE I. The nonzero improvement coefficients for 
Schrodinger functional boundary conditions with electric 
background field for various gauge groups and fermion rep- 
resentations. 

To illustrate the need for improvement, in Fig. Q] we 
show Si defined in ([8]), to compare the unimproved and 
improved results for SU(N) gauge theory with fundamen- 
tal fermions for N = 2, 3 and 4. From the figure we see 
that improvement is essential for taming the 0(a) arte- 
facts. 
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FIG. 1. Comparison between unimproved and improved re- 
sults for the fermionic contribution to the perturbative step 
scaling function for SU(2), SU(3) and SU(4) gauge theories 
with fundamental representation matter fields. 

Next, in addition to using the coefficients in Table U 
the precise values of the boundary fields must be chosen 
carefully to guarantee rapid convergence of the results. In 
Fig. [5]we show the results for S u for SU(2) gauge theory 
and adjoint fermions. The results are normalized to the 
continuum values While the improvement is necessary to 
get rid of 0(a) corrections, clearly, in order to guarantee 
that the 0(a 2 ) corrections remain small the boundary 
values must be optimized. The optimal choice is p = tt/2 
and < i] < tt/2 excluding the value r\ = 7r/4. This 
can be understood as follows: If the boundary field is 
generically a diagonal matrix with eigenvalues <f> and — 0, 
then the adjoint fermions see these eigenvalues as 2cj>, 
~2(f> and 0. In other words, since the background appears 
twice as large for adjoint fermions as for fundamental 
fermions, it seems plausible that the problems could be 
alleviated by halving the background field. We note that 
the numerical results imply that with the above choice 
of parametrizing the boundary fields, the parameter p 
plays the main role; i.e. there is a single value of p (e.g. 



7r/2 for the adjoint representation of SU(2)) for which 
the convergence is optimal and r\ can be chosen from a 
wide interval without significantly altering the result. 
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FIG. 2. Comparison between two different boundary fields 
in SU(2) gauge theory with adjoint fermions. The thin solid 
and dashed curves correspond, respectively, to improved and 
unimproved results with p — n. The thick solid (dashed) 
curves correspond to improved (unimproved) results for p — 
tt/2. 

The analysis is straightforward also for SU(3) with 
higher representation fermions. In Fig. [3] we show \5\ — 1| 
contours in the (77, p)-planc. The conventional choice 
for the fundamental representation fermions is (/?, rj) = 
(tt/3, 0). As in the case of SU(2), when translating the 
boundary matrix to the adjoint or sextet representation, 
one sees that the background field is effectively doubled. 
This again suggests trying to narrow the values of the an- 
gular parameters by a factor of two. This is indeed what 
we observe from Fig. [31 For the adjoint representation, 
shown in the upper panel, the best convergence of 5± is 
obtained for p = it/ 6 and — 4ir/9 < rj < it/9 excluding 
the value rj = —tt/Q. Similarly, as shown in the lower 
panel of the figure, for the sextet representation we find 
that the convergence is optimal for p = 677r/150 0.457T 
and -166tt/225 < -7tt/45 excluding 77 = -67tt/150. As 
for SU(2), we find that the convergence is dominantly 
controlled by the value of the p while 77 can be chosen 
from a wide interval. For the figure we used L/a = 10, 
but the results remain quantitatively similar for other 
values of L/a as well; the dependence on L/a is similar 
to the SU(2) case shown in Fig. [2] 

Finally, the improvement coefficients must be included 
consistently In Fig [4] we show comparison between the 
improved, unimproved results as well as between the re- 
sults neglecting some of the improvement coefficients. In 
particular we see that even if the Sheikholeslami-Wohlert 
coefficient c sw = 1 is implemented, but the boundary 
improvement is neglected, the results obtained under 
Schrodinger functional calculations like in [3] are pos- 
sibly still far from the proper continuum limit. Also the 
optimal values for the boundary fields used for the im- 
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FIG. 3. Convergence of |<5i — 1| as a function of p and r\ for 
adjoint (upper panel) and sextet (lower panel) fermions in 
SU(3) gauge theory. 

proved results in Fig. @]diffcr from the values used for the 
fundamental representation as we discussed earlier. We 
have used here the values of p and r\ which give optimal 
convergence for <5i as discussed in previous paragraphs. 

In this Letter we have described developments in anal- 
yses of gauge theories with higher representation matter 
fields on the lattice. Generally, we have demonstrated 
the sensitivity of the results on the correct values of the 
improvement coefficients and boundary conditions when 
using Wilson fermions. We have provided perturbative 
values of the countertcrm coefficients required for the 
analysis and emphasized the choice of the boundary con- 
ditions. Our results suggest that in studies of higher- 
representations there is need to optimize the value of the 
background field carefully This need is driven by the 
fcrmionic contribution; we have explicitly checked that 
the effect on the pure gauge contribution is negligible, on 
the level of few percents, for the cases we have considered. 
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FIG. 4. The upper panel shows the result for SU(3) with 
adjoint fermions and the lower one for SU(3) with sextet 
fermions. The thick solid and dashed curves correspond, re- 
spectively, to improved and unimproved results. The thin dot- 
ted (dash-dotted) curves correspond to the results obtained 
by setting c sw = 1 (0) and ct = (perturbative value). 
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